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We discuss the energy and vorticity spectra of turbulent superfluid '^He in the entire temperature 
range from T = 0 up to the phase transition “A point”, T\ ~ 2.17K. Contrary to classical developed 
turbulence in which there are only two typical scales, i.e. the energy injection L and the dissipation 
scales 77 , here the quantization of vorticity introduces two additional scales, the vortex core radius 
ao and the mean vortex spacing £. We present these spectra for the super- and the normal-fluid 
components in the entire range of scales from L to ao including the cross-over scale £ where the 
hydrodynamic eddy-cascade is replaced by the cascade of Kelvin waves on individual vortices. At 
this scale a bottleneck accumulation of the energy was found earlier at T = 0. We show that even 
very small mutual friction dramatically suppresses the bottleneck effect due to the dissipation of the 
Kelvin waves. Using our results for the spectra we estimate the Vinen “effective viscosity” v' in the 
entire temperature range and show agreement with numerous experimental observation for v'{T). 


I. INTRODUCTION 

Superfluidity was discovered by Kapitza and by Allen 
and Misener in 1938 who demonstrated the existence of 
an inviscid fluid flow of ^He below Tx ~ 2.17K. In the 
same year, London linked the properties of the superfluid 
"^He to Bose-Einstein condensation. 

Soon after, Landau and Tisza offered a “two fluid” 
model in which the dynamics of the superfluid ^He is 
described in terms of a viscous normal component and an 
inviscid superfluid component, each with its own density 
Pn{T) and Ps{T) and its own velocity field u^{r,t) and 
Us{r,t). Already in 1955, Feynman realized— that the 
potential appearance of quantized vortex lines will result 
in a new type of turbulence, the turbulence of superfluids. 
The experimental verification of this prediction followed 
in the paper by Hall and Vinen a year later—. 

An isolated vortex line is a stable topological defect 
in which the superfluid density drops to zero and the 
velocity Vg oc 1/r. Here r > oq is the radial distance 
from the center that exceeds a core radius ao — 10“® cm 
in ^He. The existence of quantized vortex lines^^— in 
superfluid tubulence introduces automatically additional 
length scales that do not exist in classical turbulence. In 
addition to Oq, the density of vortex lines £ defines an 
“inter-vortex” average spacing denoted as £ = 1/ x/£. 

The pioneering experimental observation of Maurer 
and Tabeling— showed quite clearly that the large-scale 
energy spectrum of turbulent "‘He above and below Tx 
are indistinguishable. This and later experiments and 
simulations gave rise to the growing consensus that on 
scales much larger than I the energy spectra of turbu¬ 
lent superfluids are very close to those of classical fluids 
if they are excited similarly—. The understanding is that 
motions with scales £ correlate the vortex lines, or¬ 
ganizing them into vortex bundles. At these large scales 
the quantization of the vortex lines becomes irrelevant 
and the large scale motions are similar to those of con¬ 
tinuous hydrodynamic eddies. Obviously, since energy is 
cascaded by hydrodynamic eddies to smaller and smaller 


scales, we must reach a scale where the absence of viscous 
dissipation will require new physics. 

Evidently, when the observation scales approach i and 
below, the discreteness of the quantized vortex lines be¬ 
comes crucial. Indeed, on such scales the dynamics of the 
vortex lines themselves become relevant including vortex 
reconnections and the excitation of Kelvin waves on the 
individual vortex lines. Kelvin waves exist also in clas¬ 
sical hydrodynamics but here they become important in 
taking over the role of transferring the energy further 
down the scales. Their nonlinear interaction results in a 
so-called “weak wave turbulence” — 1 ^ supporting a mean 
energy flux towards shorter and shorter scales. Finally 
when the cascade reaches the core radius scale the energy 
is radiated away by quasiparticles (phonons in ^He )— . 

Although the overall picture of superfluid turbulence 
described above seems quite reasonable, some impor¬ 
tant details are yet to be established. A particularly 
interesting issue is the physics on scales close to the 
crossover between the eddy-dominated and the Kelvin 
wave-dominated regimes of the spectrum. It was pointed 
out in Ref.— that the nonlinear transfer mechanism of 
weakly nonlinear Kelvin waves on sparse vortex lines is 
less efficient than the energy transfer mechanism due to 
strongly nonlinear eddy interactions in continuous flu¬ 
ids. This may cause an energy cascade stagnation at the 
crossover scale. 

The present paper is motivated by some exciting new 
experimental and simulational developmentsi^^— that 
call for a fresh analysis of the physics of superfluid turbu¬ 
lence in a range of temperatures and length scales. These 
developments include, among others, cryogenic flow visu¬ 
alization techniques by micron-sized solid particles and 
metastable helium molecules, that allow, e.g. direct ob¬ 
servation of vortex reconnections, mean normal and su¬ 
perfluid velocity profiles in thermal counterflow — the 
observation of Andreev reflection by an array of vibrating 
wire sensors shedding light on the role of vortex dynam¬ 
ics in the formation of quantum turbulence, etc.—; the 
measurements of the vortex line density by the atten¬ 
uation of second soundi^^— and by the attenuation of 
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ion beam a^^i^^ . An important role in the recent develop¬ 
ments is played by large-scale well-controlled apparata, 
like the Prague^^^— and the Grenoble wind tunnels^^ , 
the Manchester spin-dow n^^i^° , the Grenoble Von Kar- 
man flowa^i, the Helsinki rotating cryostat^^— , and 
some other experiments. Additional insight was provided 
by large-scale numerical simulations of quantum turbu¬ 
lence by the vortex-filament and other methods that gives 
direct access to detailed picture of vortex dynamics which 
is still unavailable in experiments, see also Refs.— . 

The stagnation of the energy cascade at the intervor¬ 
tex scale mentioned above is referred to as the bottle¬ 
neck effect. This issue was studied in Ref.— in the ap¬ 
proximation of a “sharp” crossover. Ref.— introduced 
a model of a gradual eddy-wave crossover in which both 
the eddy and the wave contributions to the energy spec¬ 
trum of superfluid turbulence at zero temperature £^{k) 
(see Eq. (llOaD l were found as a continuous function of 
the wave vector k. The main message of Ref.— is that 
the bottleneck phenomenon is robust and common to all 
the situations where the energy cascade experiences a 
continuous-to-discrete transition. The details of the par¬ 
ticular mechanism of this transition are secondary. In¬ 
deed, most discrete physical processes are less efficient 
than their continuous counterparts^!. On the other hand, 
particular mechanisms of the continuous-to-discrete tran¬ 
sition can obviously lead to different strengths of the bot¬ 
tleneck effect. 

The main goal of the present paper is to develop a 
theory of superfluid turbulence that analyzes the dy¬ 
namics of turbulent superfluid ^He and computes its 
energy and vorticity spectra in the entire temperature 
range from T —^ 0 up to the phase transition, “A point” 
Tx ~ 2.17 K, and in the entire range of scales R, from 
the outer (energy-containing, or energy-injection) scale 
L down to the core radius oq. We put a particular fo¬ 
cus on the crossover scales R ^ i, where the bottleneck 
energy accumulation is expected— 

The main results of this paper are presented in Sec. m 
Its introductory subsection, 

III Al Basic approximations and models, 

overviews the basic physical mechanisms, which 
determine the behavior of superfluid turbulence 
and describes the set of main approximations and 
models, used for their description. 

The rest of Sec. m is devoted to the following 
problems: 

III B1 Temperature dependence of the energy spectra and 
the bottleneck effect in turbulent ^He; 

III Gl Temperature dependence of the vorticity spectra; 

nmi Gorrelations between normal and superfluid mo¬ 
tions and the energy exchange between compo¬ 
nents; 


III El Temperature dependence of the effective superfluid 
viscosity in ^He. 

Glearly, the basic physics of the large scale motions, 
differ from that of small scale motions. The same can 
be said about different regions of temperature: zero tem¬ 
perature limit, small, intermediate and large tempera¬ 
tures. It would be difficult to follow the full description of 
the physical picture of superfluid turbulence in all these 
regimes without clear understanding of the entire phe¬ 
nomenon as a whole. Therefore in Sec. lII Al we restricted 
ourselves to a panoramic overview of the main approx¬ 
imations and models, leaving detailed consideration of 
some important, but in some sense secondary issues, to 
the next two sections of the paper fSec. lIIIl and Sec. lIVI) . 
These include the analysis of the range of validity of the 
basic equations of motions, of the main approximations 
made in the derivation, and of the numerical procedures. 
These sections consist of the following subsections: 

IIII Al Goarse-grained, two-fluid, gradually-damped Hall- 
Vinen-Bekarevich-Khalatnikov (HVBK) equations; 

IIII Bl Two-fluid Sabra shell-model of turbulent ^He; 

IIV Al Differential approximations for the energy fluxes 
of the hydrodynamic and Kelvin wave motions; 

IIVBI One-fluid differential model of the graduate eddy- 
wave crossover; 

In the final Section El we summarize our results on the 
temperature dependence of the energy spectra of the 
normal and superfluid components in the entire region 
of scales. We demonstrate in Eig.[5]that the computed 
temperature dependence of the effective viscosity v'iT) 
agrees qualitatively with the experimental data in the en¬ 
tire temperature range. We consider this agreement as a 
strong evidence that our low-temperature, one fluid dif¬ 
ferential model and the high temperature coarse-grained 
gradually damped HVBK model capture the relevant ba¬ 
sic physics of the turbulent behavior of ^He. 

II. UNDERLAYING PHYSICS AND THE 
RESULTS 

A. Basic approximations and models 

1. Coarse-grained, two-fluid, gradually-damped HVBK 
eguations 

As we noticed in the Introduction, the large-scale mo¬ 
tions of superfluid ^He (with characteristic scales R:^ £) 
are described using the two-fluid model as interpene¬ 
trating motions of a normal and a superfluid compo¬ 
nent with densities pn{r,t,T), ps{r,t,T) and velocities 
Un{r,t), Us{r,t). Following Ref.— we neglect variations 
of densities by considering them as functions of the tem¬ 
perature T only, Pn{T) and Ps{T). We also neglect both 
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the bulk viscosity and the thermal conductivity. This re¬ 
sults in the simplest form of the two incompressible-fluids 
model for superfluid ^He that have a form of the Euler 
equation for itg and the Navier-Stokes equation for 
see e.g. Eqs. (2.2) and (2.3) in Donnely’s textbook^. As 
motivated below, we add an effective superfluid viscosity 
term also in the superfluid equation, writing 

d u 1 

+ {Us-V)Us -Vps = i/'Alts-Fns, (la) 

Ot Ps 

d u 1 

-k (itn ■ V)Mn-Vpn = t'n AWn + —i^ns ■ (lb) 

ot Pn Pn 

Here Pn, Ps are the pressures of the normal and the su¬ 
perfluid components: 

Pn = ^[p + Ps\Us - Un\'^] , Ps = ^ [p - - Mn^] , 

p = Ps + Pn is the total density, z/n is the kinematic vis¬ 
cosity of normal fluid. 

The term describes the mutual friction between 
the superfluid and the normal components mediated by 
quantized vortices, which transfer momentum from the 
superfluid to the normal subsystem and vice versa. Fol¬ 
lowing Ref.—, we approximate it as follows: 

Fug — ot PgCdgi^Ug "Wn) ; 

where uig is the characteristic superfluid vorticity. 

The equations 0 are referred to as the Hall-Vinen- 
Bekarevich-Khalatnikov (or HVBK) coarse-grained 
model. The relevant parameters in these equations, 
are the densities Ps{T) and Pn{T), the mutual friction 
parameters ot[T) and the kinematic viscosity of the 
normal-fluid component Vn{T) normalized by pn. 

The original HVBK model does not take into account 
the important process of vortex reconnection. In fact, 
vortex reconnections are responsible for the dissipation 
of the superfluid motion due to mutual friction. This ex¬ 
tra dissipation can be modeled as an effective superfluid 
viscosity VgiT) as suggested in Ref.—: 

Vg{T)^aK. (Id) 

We have added a dissipative term proportional to v'g to 
the standard HVBK model and the resulting Eqs. ([T|) [dis¬ 
cussed in more details in Sec. m will be referred to as 
the “gradually damped HVBK model”. We use this name 
to distinguish our model from the alternative “truncated 
HVBK model” suggested in Ref.— which was recently 
used for for numerical analysis of the effective viscosity 
it'{T) in Ref.—. We suspect that the sharp truncation 
introduced in the latter model creates an artificial bot¬ 
tleneck effect that is removed in the gradually damped 
model. The difference in predictions between the models 
will be further discussed in Sec. lHI A'Sl 

2. Two-fluid Sabra shell-model of turbulent *He 

The gradually damped HVBK Eqs. © provide an ad¬ 
equate basis for our studies of the large-scale statistics 


of superfluid turbulence. However their mathematical 
analysis is very difficult because of the same reasons that 
make the the Navier-Stokes equations^ difficult. The in¬ 
teraction term is much larger than the linear part of the 
equation (their ratio is the Reynolds number, Re^ 1), 
the nonlinear term is nonlocal both in the physical and 
in the wave-vector fc-space, the energy exchange between 
eddies of similar scales, that determines the statistics of 
turbulence, is masked by much larger kinematic effect of 
sweeping of small eddies by larger ones, etc. 

Direct numerical simulations of the HVBK Eqs. © are 
even more difficult than the Navier-Stokes analog, being 
extremely demanding computationally, allowing there¬ 
fore for a very short span of scales. A possible simpli¬ 
fication is provided by shell models of turbulence— “—. 
They significantly simplify the Navier-Stokes equations 
for space-homogeneous, isotropic turbulence of incom¬ 
pressible fluid. The idea is to consider the equations in 
wave vector fc-Fourier representation and to mimic the 
statistics of u{k,t) in the entire shell of wave numbers 
km < km+i by only one complex shell velocity Vm- The 
integer index m is referred to as the shell index, and the 
shell wave numbers are chosen as a geometric progression 
km = koX^, with A being the shell-spacing parameter. 
This results in the ordinary differential equation 

(s = NLm{?;m4 . (2a) 

Here the nonlinear term NLm{r'm} is linear in k and 
quadratic in Vm' ( a functional of the set {v!^}), which 
usually involves shell velocities with |m — m'| < 2 . the 
kinetic energy is preserved by the nonlinear term. For ex¬ 
ample, in the popular Gledzer-Ohkitani-Yamada (GOY) 
shell model^i^ 

= zA:^(nAu^_|_2U^_|_i -\-bv^—iV^^i 
+CVm-lVm-2r , GOY , (2b) 

where the asterisk * stands for complex conjugation. In 
the limit —?► 0 and with a -|- 6 -I- c = 0, Eqs. ([2|) preserve 
the kinetic energy E = and has a second in¬ 
tegral of motion H = c)'^\vm\^. The traditional 

choice a = A|6| allows to associate H with the helicity in 
the Navier-Stokes equations. 

Note that the simultaneous rescaling a => ap, b ^ bp 
and c ^ cp with some factor p results in a straightfor¬ 
ward rescaling of the the time variable t ^ t/p without 
any effect on the instantaneous stationary statistics of 
the model. Thus, the shell model ([2]) has only one fitting 
parameter A, which has only little effect on the resulting 
statistics. The traditional choice A = 2 allows to reason¬ 
ably model the interactions in fc-space with an efficient 
energy exchange between modes of similar index m. 

We stress that with the above choice of parameters, 
a-|-6-|-c = 0, a = A|6|, and 

A = 2, 0 = 1, b = c = —0.5 , 


(3) 
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the shell models reproduced well various scaling prop¬ 
erties of space-homogeneous, isotropic turbulence of in¬ 
compressible fluids, see Ref.— and references therein. To 
mention just a few: 

- the values of anomalous scaling exponents (see, e.g. 
Table I ini^); 

- the viscous corrections to the scaling exponents—; 

- the connection between extreme events (outliers) and 
multiscaling—; 

- the inverse cascade in two-dimensional turbulence—; 

- the strong universality in forced and decaying 
turbulence;^, etc—“— . 

Therefore, we propose shell models are a possible alter¬ 
native to the numerical solution of the HVBK Eqs. ([T|). 
This option was studied in Ref.— which proposed a two- 
fluid GOY shell model for superfluid turbulence with an 
additional coupling by the mutual friction. 

In our studies of superfluid turbulence— and below, 
we use the so-called Sabra-shell model—, with a different 
form of the nonlinear term: 


-cVm-iVm-2 ), Sabra . (4) 

The advantage of the Sabra model over the GOY model 
is that the resulting spectra do not suffer from the un- 
physical period-three oscillations, thanks to the strong 
statistical locality induced by the phase invariance— 

We solved numerically the two-fluid Sabra-shell model 
form of the HVBK equations (l2a)) and (H)) coupled by the 
mutual friction, for the shell velocities. Gathering enough 
statistics, we computed the pair- and cross-correlation 
functions of the normal- and the super-fluid shell veloc¬ 
ities. This led to the energy spectra 6^ {k) and £s (k) 
together with the cross-correlation £ns{k). 

In the simulations we used 32 shells. All the results are 
obtained by averaging over about 500 large eddy turnover 
times. The rest of details of the numerical implementa¬ 
tion and simulations are given in Sec. HUB] 


3. Low temperature one-fluid eddy-wave model of 
superfluid turbulence 

As we just explained, in the high-temperature region 
the fluid motions with scales R < £ are damped and 
motions with R> £ are faithfully described by the Sabra- 
shell model (|2a1) and (H]) . In this approach we first solve 
the dynamical equation and then perform the statistical 
averaging numerically. 

In the low temperature regime, T < T\/2, where the 
Kelvin wave motions of individual vortex lines are impor¬ 
tant this approach is no longer tenable. Instead, we adopt 
a different strategy, in which we first perform the statis¬ 
tical averaging analytically and then solve the resulting 
equations for the averaged quantities numerically. 

To this end we begin with the dynamical Biot-Savart 
equation of motion for quantized vortex lines. Then we 


applied the Hamiltonian description to develop a “weak 
turbulence” formalism to the energy cascade by Kelvin 
waves^. This approach results in a closed form expression 
for the Kelvin wave energy spectra, derived in Ref. 

^ ’ LN-spectrum. (5a) 

Here Ekw is the energy flux over small-scale region, R < £, 
and A ~ ln(£/ao). The value of the universal constant 
Cln ~ 0.304 was estimated analytically in Ref. HI. The 
dimensionless constant 4/ may be considered as the r.m.s. 
vortex line deflection angle at scale £ and is given by 


. T , _ £kw£'^ 

= Ak2 


(5b) 


In the low-temperature region, T < Tc/2, the den¬ 
sity of the normal component is very small and due to 
very large kinematic viscosity it may be considered at 
rest. Therefore the large scale motions of ^He, R > £, 
are governed by the first of HVBK Eq. (llaL which coin¬ 
cides with the Navier-Stokes equation in the limit T —>■ 0. 
Therefore, in the hydrodynamic range of scales, R > £, 
we can use the Kolmogorov-Obukhov 5/3-law^ for the 
hydrodynamic energy spectrum: 

fHD(fc) = , KO-41 . (6) 


Here Ehd is the energy flux over large scale range and 
Ckh ~ 1 is the Kolmogorov dimensionless constant. 

Both spectra, ([5]) and ([5]) have the same /c-dependence, 
oc but different powers of the energy flux. A way 

to match these spectra in the T —>• 0 limit was suggested 
in Ref.—. The idea was to adopt the differential approxi¬ 
mations to the Kelvin-wave^ and hydrodynamic-energy 
flux—, based on their spectra (153 and (O : 


S^HD {k) 


£kw(^) 


, Pa) 

3f2^(A:)vI/2A:6£4 d£^^{k) 

5(ClnAk)3 dk ’ ^ ^ 


and to construct a differential approximation for the su¬ 
perfluid energy flux e that is valid for all wave numbers 
(including the cross over scale): 

KW HD 

e(fc) = eHD(fc) + eKw(fc) + Ehd (^) + £kw(^) • (^c) 

KW HD 

The additional cross-contributions (fc) and Sy^^{k) 

originate from the interaction of two types of motion, 
hydrodynamic and Kelvin waves. 

Eor r —0 the total energy flux should be k- 
independent, e{k) =const. As explained in Ref.— this 
leads to an ordinary differential equation for the total 
superfluid energy £s{k) = ^’hd)^) + £K.w(k). In this pa¬ 
per we generalize this approach to the full temperature 
interval with the help of the energy balance equation 


d£s{k,t) 

dt 


des{k,t) 

dk 


= v'^k'^£s{k,t) . 


( 8 ) 
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The right-hand-side of this equation originates from the 
Vinen-Niemella viscosity in Eq. dial) and accounts for the 
dissipation in the system. 

Much more detailed description of this procedure can 
be found in Sec.lIVl 

B. Temperature dependence of the energy spectra 
and the bottleneck effect in turbulent ^He 

To discuss our results we define the energy spectrum 
£(k) of isotropic turbulence (in one-dimensional fc-space) 
such that 

OO 

^ = j ^{k)dk=]^{\u{r,t)\^) , (9) 

0 

is the energy density £ of ^He per unit mass. Hereafter 
(...) stands for the “proper averaging” which may be 
time averaging over long stationary dynamical trajectory 
or/and space averaging in the space-homogeneous case, 
or an ensemble averaging in the theoretical analysis. As¬ 
suming that the Navier-Stokes dynamics are ergodic, all 
these types of averaging are equivalent. 

In the low temperature range, where ^He consists 
mainly of the superfluid component, we distinguish the 
spectrum of large scale hydrodynamic motions with k(. 

1, denoted as £^^{k), from the spectrum of small scales 
Kelvin waves (with k^ ^ 1), denoted as £^^{k). The 
total superfluid energy spectrum is written as 

= . (10a) 

In the high temperature range, where the densities of 
the super-fluid and normal-fluid components are compa¬ 
rable, but Kelvin waves are fully damped, we will distin¬ 
guish the spectrum of hydrodynamic motions of the su¬ 
perfluid component at large scales as f®(fc), from that of 
the normal-fluid component, f"(/c), omitting for brevity 
the subscript “hd” • In this temperature range, the total 
energy spectrum of superfluid ^He is written as 

£{k) = £^{k)+£'^{k) . (10b) 

The resulting energy spectra ^’^□(fc), £^{k) and £’^{k) 
for a set of eleven temperatures from T = 0.32 K to 
T = 2.16 K are shown in Fig.[T] 


1. Low-temperature one-fluid energy spectra 

First, we discuss the results for the eddy-wave model 
of superfluid turbulence [cf. Sec. Ill A 3j . for the low 
temperature range T < T\/2 ~ 1.08 K, which is shown 
in Fig.[T^. These spectra are compensated by a fac¬ 
tor {ktff^ such that both the Kolmogorov-Obukhov- 
41 spectrum £^j^{k) oc Eq. dH]) (for the hydro- 

dynamic scales kt 1) and the Lvov-Nazarenko spec¬ 
trum £^^{k) oc Eos. (153 show up as a plateau. 


These plateaus are clearly seen for the lowest shown tem¬ 
perature T = 0.32 K. Moreover, the full energy spec¬ 
trum (solid blue line) £^{k) demonstrates the existence 
of an important bottleneck energy accumulation. We ob¬ 
serve a large cross over region connecting the HD re¬ 
gion /t£ < 1, where £^{k) —>■ £^^(k) with a much higher 
plateau £^„(k) for k£ > 50 where £‘^(k) —>■ £^^(k). 

In the cross over region the compensated energy spec¬ 
trum is close to (fcl')®/^ (cf. the black dashed line), mean¬ 
ing that £^(k) depends on k only weakly. In this re¬ 
gion the energy spectrum is dominated by Kelvin waves, 
£^(k) ~ £l®j,(fc), while the energy flux in dominated by 
the HD eddy motions. Therefore we have here a flux-less 
regime of Kelvin waves. Without flux the situation re¬ 
sembles thermodynamic equilibrium, in which the Kelvin 
waves energy spectrum corresponds to energy equiparti- 
tion between the degrees of freedom, i.e. iJ®^(fc) ^const, 
as observed. 

For k£ > 2 ■ 10^ the Kelvin waves energy spectrum 
at T = 0.32 K is suppressed by the mutual friction, as 
explained in Sec. lIV Al In FiglT^ this part of the spectrum 
is not shown; however one sees progressive suppression 
of the energy spectra with temperature increasing from 
0.44 K (around k(. ~ 400) to T = 1.07 K (around ~ 2). 
It is important to notice that for T > 0.49 K the HD part 
of the spectrum is practically temperature independent; 
only the Kelvin waves energy spectra are suppressed by 
the temperature, cf. the coinciding dashed lines in Fig.[T]i 
for T = 0.32 , 0.44 and 0.49 K. 

For T > 0.5 K, the Kelvin wave contributions to the 
energy spectra are very small - the solid and the dashed 
lines for the same temperature are fairly close. Finally, 
the dashed and the solid lines for T = 1.07 K practically 
coincide, i.e. the Kelvin waves are fully damped. This 
means that for T > 1 K there is no need to account 

for the Kelvin wave motions on individual vortex lines, 
and the full description of the problem is captured by the 
coarse-grained HVBK. 


2. High-temperature two-fluid energy spectra 

The energy spectra, obtained with the Sabra-shell 
model form of HVBK equations (Eal) and for tem¬ 
peratures T > T\/2 ~ 1.08 K, are shown Fig.jTjD for 
T = 1.07K (in blue), T = 1.3K (in magenta), 1.8K (in 
green) and T = 2.16 K (in red). The lowest temperature 
in this two-fluid approach, T = 1.07 K, was chosen for 
comparison with the highest temperature 1.07 K in the 
one-fluid approach; see Fig.[T]a,. At T = 1.3 K, which is a 
frequently used temperature in numerical simulations of 
superfluid turbulence, the normal fluid component is not 
negligible (pn/p — 0.045), and the normal fluid kinematic 
viscosity is still much larger than that of the superfluid: 
v^/v'^ ~ 23. For T = 1.8, when pi^jp ~ 0.3, the kinematic 
viscosities are close to each other (see Tab.lJ). At higher 
temperatures the normal fluid components play more and 
more important role until they dominate at T > 2.0 K, 
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T, K 

0.43 

0.55 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2.0 

2.1 

2.16 

Pn/p 

- 

- 

9.3 • IQ-'^ 

0.003 

0.007 

0.014 

0.026 

0.045 

0.0728 

0.111 

0.162 

0.229 

0.313 

0.420 

0.553 

0.741 

0.907 

a 

10 "® 

10 “® 

6.5 • lO”"* 

0.0025 

0.0056 

0.011 

0.026 

0.034 

0.051 

0.072 

0.097 

0.126 

0.160 

0.206 

0.279 

0.48 

1.097 

ap/p^ 

- 

- 

0.70 

0.83 

0.80 

0.78 

1.00 

0.76 

0.70 

0.65 

0.60 

0.55 

0.51 

0.49 

0.53 

0.65 

1.209 

V ! K 

- 

- 

1.09 

0.43 

0.27 

0.17 

0.12 

0.10 

0.10 

0.09 

0.09 

0.09 

0.09 

0.093 

0.101 

0.124 

0.154 

ZZn/zt 

- 

- 

1179 

148 

38 

11.1 

4.62 

2.34 

1.32 

0.84 

0.56 

0.39 

0.29 

0.22 

0.182 

0.167 

0.170 

v'Jk 

- 

- 

0.0067 

0.022 

0.040 

0.061 

0.099 

0.101 

0.135 

0.171 

0.207 

0.234 

0.237 

0.280 

0.312 

0.427 

0.815 


TABLE I: The parameters of the superfluid ^He, taken from Refs.— the relative density of the normal component Pn/p, 
the mntual friction parameter a, the combination apjpn [which weakly depends temperature and is responsible for the mutual 
friction density in Eq. O], He-II kinematic viscosity v = p/p {p is the dynamic viscosity) and the kinematic viscosity of the 
normal-fluid component i/n = p/pn', the effective superfluid viscosity (inter- and extrapolation of Ref.— results). 

(a) Low T < Tx/2, one-fluid model (b) High T > Ty/2, two fluid model 




FIG. 1: Color online. Log-Log plots of the energy spectra in superfluid ^He compensated by the inertial range scaling at 
different temperatures (shown as labels). Panel (a): Plots of the (compensated by and normalized by their values at 
kl — 0.1) full superfluid energy spectra (solid lines) and their hydrodynamic (large scale) parts £hd(«^^) (dashed lines) 

for the one-fluid model fSec. lIVB]| . Panel (b): Plots of the (compensated by the anomalous scaling and normalized by 

their inertial range value) shell energies of the normal fluid component = k£^(km£) (solid lines) and of the superfluid 

component = k£’^{kmi) (dashed lines) for the two-fluid shell model fSec. lIII Bf . 


when pn > Ps- At the highest temperature in this simula¬ 
tion, T = 2.16 iC, close to T\, we have pn ~ 0.9p, and the 
effective superffuid kinematic viscosity z/' is even larger 
than i/n. 

Shell-model simulations reproduce intermittency ef¬ 
fects and therefore the scaling exponent ^2 of the en¬ 
ergy spectra £{k) oc slightly differs from the KO- 
41 prediction, ^2 7 ^ 5/3. For the chosen shell-model 
parameters^Siffi ^2 ~ 1-72 which is quite close to the ex¬ 
perimental observations. For better comparison with the 
low-temperature one-fluid results of Fig.[T^, we show in 
Fig.Hb the normal (solid lines) and superfluid (dashed 
lines) energy spectra £''{k£) and £^{kl), compensated by 
(M)^'^^ so that they exhibit a plateau in the inertial in¬ 
terval of scales. 


As expected, for T = 1.8 K, when ~ r's, the super¬ 
fluid and normal fluid spectra are very close, and similar 
to the spectra of classical fluids. In the inertial range they 
demonstrate the anomalous behavior £1® oc oc 
[|WmP oc k^ with the scaling exponent .^2 ~ 1.72. 
Moreover, due to the strong coupling between the normal 
and superfluid component (discussed below in Sec. Ill D() 
the energy fluxes in both components are equal (see, e.g. 
Fig. HI), and therefore the energies are equal in the iner¬ 
tial interval as well, £’®(fc) = £^{k). Non-trivial behavior 
occurs only in the inertial-viscous crossover region; there¬ 
fore the inertial interval is not shown in Fig.[TJ 

For T = 1.07, when ~ 180 zz', the viscous cutoff 
of the normal fluid’s spectrum, /cmax.n, occurs at much 
smaller k than the cutoff of the superfluid spectrum. 
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(a) Low T < Ta/2, one-fluid model 


(b) High T > Ta/ 2, two fluid model 




FIG. 2: Color online. Color online. Linear-log plots of the vorticity spectral densities, {k£) Q,^{k£), normalized by their k£ = 0.1 
values, at different temperatures. The total mean-square vorticity (w^) is proportional to the area under the plot. Panel (a): 
Superfluid vorticity spectra {k£) Q,^{k£) at low temperatures in the one-fluid model, Sec. lIV Bl Panel (b): Normal and superfluid 
vorticity spectra {k£) Q^{k£) (solid lines) and {k£) ^^{kl) (dashed lines) at high temperatures in the two-fluid model, Sec. lIII Bl 


fcinax.s- To estimate the ratio fcmax.s/fcmax.n, notice that 
in the KO-41 picture of turbulence fcmax may be found 
by balancing the eddy-turnover frequency, 

7(fc) ~ eT3;,2/3 _ ^3/2^^^ 

with the viscous dissipation frequency This gives 

the well known result 

fcmax - . (12) 

In our case Eg = En- Therefore, neglecting the energy 
exchange between the super- and the normal-fluid com¬ 
ponents, we get an estimate: 

^max,n ^ ‘^s '' 

For T = 1.07, when — 180 zz' this gives 

fcmax.s/fcmax,n — 50—in a good agreement with the re¬ 
sult in Fig. [lb. For T = 1.3 K, the ratio of the viscosities 
is smaller (about ~ 23, see Tab. HI . Therefore the dif¬ 
ference in cutoffs is less pronounced. As expected, for 
T = 2.16 K, when zz^ ~ 0.2zz' the situation is the oppo¬ 
site, and the superfluid component is damped at a smaller 
k than the normal one. 

Notice that there is no bottleneck energy accumula¬ 
tion in the spectra (see Figs. jib) obtained using the shell 
model approximation of the gradually damped HVBK 
equations. This is qualitatively different from the results 
of the truncated HVBK model—, which demonstrated a 
very pronounced bottleneck both in the normal and the 
superfluid components, e.g. at T = 1.15 K. The latter 
would lead to a huge contribution to the mean square 
superfluid vorticity (jwsp) and, as a result, to a very 


small effective Vinen’s viscosity zz'. This would definitely 
contradict the experimental observation shown in Fig.jS] 
We will discuss this issue in greater detail in Sec. Ill El 

C. Temperature dependence of the vorticity 
spectra in turbulent “^He 

At this point we cannot compare our predictions for 
energy spectra with experimental observations, especially 
in the cross-over and in the small scale regions. This 
stems from the lack of small probes, see cf. the review^. 
On the other hand, the attenuation of second sound or 
ion scattering may be used to measure the mean vortex 
line density £ = in '^He or even its time and space 
dependence^. In turn, the value can be expressed in 
terms of the mean-square superfluid vorticity (lazgp) via 
the quantum of circulation 

(|a;g|2) « {nCf . (14) 

Therefore, the information about the vorticity is very im¬ 
portant from the viewpoint of comparison with available 
and future experiments. 

By analogy with the energy spectra (lH, let us define 
the power spectra of vorticity n^(fc) so that the mean- 
square vorticity (|wp) is given by the integral: 

CO CO 

(|w|2) = Jn^{k)dk = Jk n^{k)d\nk . (15) 

0 0 

In isotropic incompressible turbulence = 

2k'^£-au{k)- Therefore we define 

nl{k) = 2k'^£l^{k ), nl{k) = 2p£’^{k). (16) 





























For brevity, we omit the subscript “hd” for the normal 
component; 

Plots of kVilj^{k) for different temperatures are shown 
in Figured] According to Eq. (ITSl) . the area under these 
plots is proportional to the total mean square vorticity, 
(|ws,nP)- Fig.d^ shows the results for the eddy-wave 
model (the corresponding energy spectra for the same 
temperatures are shown in Fig.[T]a,). One sees that the 
largest (and temperature independent) value of (jccsp) 
is reached for T < 0.49 K: plots for T = 0.32 0.44 and 
0.49 K practically coincide. Accordingly, the temperature 
range T < 0.49 K may be considered as zero-temperature 
limit with the maximal value of (and correspond¬ 

ingly, the smallest value of r's(F), as we will discuss later). 
At temperatures above 0.5K the area under the plots de¬ 
creases (and correspondingly, Vs{T) increases). 

In Fig.db we show vorticity spectra {ktj (|a;s,n ?m) 
of the normal-fluid (solid lines) and the superfluid compo¬ 
nents (dashed lines) for different temperatures obtained 
in the framework of the Sabra-shell model (the corre¬ 
sponding spectra are shown in Fig.[Tj3). Again, the area 
under the plots is proportional to the total mean square 
vorticity (Iwn.sP)- One sees that for the lowest temper¬ 
ature T = 1.07 K the normal fluid vorticity (blue solid 
line) is fully suppressed by the huge normal viscosity, 
while the superffuid vorticity is very large. At this tem¬ 
perature one can describe the superfluid "‘He in the range 
of scales ~ 1 using a one-fluid approximation with zero 
normal-fluid velocity. This provides the main contribu¬ 
tion to the vorticity. To some extent, this situation per¬ 
sists up to T Ri 1.3 K, when the superffuid vorticity is still 
larger than the normal one, see Fig.[T]3. As expected, for 
T ~ 1.8, when the normal and superffuid viscosities are 
compatible, the normal and superffuid vorticities are very 
close. For these and higher temperatures the analysis of 
our problems definitely calls for a two-fluid description. 


D. Correlations of normal and superfluid motions 
and energy exchange between components 

1. Correlations of the normal and superfluid velocities 

It is often assumed (see e.g. RefJ^) that the normal 
and superffuid velocities are “locked” in the sense that 

Un{r,t) = Us{r,t) , (17) 

(at least in the inertial interval of scales). For quan¬ 
titative understanding to which extent this assumption 
is statistically valid we consider the simplest possible 
case of stationary, isotropic and homogeneous turbulence. 
Here we introduce a cross-correlation function (in ID 
fc-representation) of the normal and the superffuid ve¬ 
locities fns(fc)- This correlation function is defined us¬ 
ing the simultaneous, one-point cross-velocity correlation 
{un{r,t) ■ Us{r,t)) similarly to Eq. ([5]): 

J £ns{k)dk = ^{un{r,t)-Us{r,t)) . (18) 


High T > Tx/2, two-fluid model 



FIG. 3: Color online. Cross-correlation coefficients 
Eq. (I20al) (solid lines) and K 2 {ki) , Eq. (I20bll (dashed lines) 
for different temperatures. Color code is the same as in Fig.[T]D 
and Fig.[2h: T = 1.07 K - blue, T — 1.3 K - cyan, T = 1.8 K 
- green, and T = 2.16 K - red. 


If, for example, motions of the normal and the superffuid 
components at a given k are completely correlated, then 
£ns{k) = £n{k) = £s{k). If this is true for all scales, then 
Eq. (ITTl) is valid. 

It is natural to normalize fns by the normal and the 
superffuid energy densities, £n and £s. This can be rea¬ 
sonably done in one of two ways: 


or 


JCi{k) 
IC2{k) E 


2£'ns(fc) 

^ £n{k)+£,{k) ’ 

£ns{k) 

x/£n{k) - £,{k) ■ 


(I9a) 

(19b) 


Both coefficients are equal to unity for fully locked su¬ 
perffuid and normal velocities, Eq. (II3, and both van¬ 
ish if the velocities are statistically independent. How¬ 
ever, if Un{r,t) = Cus{r,t), with C ^ 1 then /Ci(fc) = 
2CI{C'^ 4-1) < 1, but still JC 2 {k) = 1. In any case 
/Ci(fc) ^/C2(fc). 

In shell models, the coefficients /Ci(fc) and K. 2 (k) can 
be written as follows: 


Ah(fc™£) 

K2{km^) 


2Re(<*<) 

’ 


(20a) 

(20b) 


These objects are shown in Eig.[3| At first glance, it is 
surprising that the correlations K 2 {k^f) (dashed lines in 
Fig-Ej) for T ^ 1.8 K persist for much larger wave vectors 
than Ki{km(-), approaching km.(- ~ 10^. For example, for 
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T = 1.07K (blue lines) Ki{km() vanishes at ~ 1, 
while K 2 {kmtj > 0.95 all the way up to kmi — 100. In 
this range of scales (1 < km-^ < 100) <C but 
v^{t) oc u^(t), meaning that strongly damped normal 
velocity does not have its own dynamics and should be 
considered as “slaved” by the superfluid velocity. The 
damped velocity (normal or superfluid) at any tempera¬ 
ture T ^ 1.8K would follow this ’’slaved” dynamics. 

A model expression of the cross-correlation fgn in 
terms of the self-correlation functions Eg and £n was 
found in Ref.—. In current notations it reads: 


fsn(fc) 


_ aUJs[pn£n{k) -I- PsEsjk)] _ 

aUJsP + Pn[(r's + l^n) + 7 n(fc) + 7 s (fc)] 


,( 21 ) 


where the characteristic interaction frequencies (or 
turnover frequencies) of eddies in the normal and super¬ 
fluid components, 7n(^) and 7s (fc), are given by Eq. (fTTl) 
and Qg is defined as: 

(7, = V(|wsP) . (22) 


The derivation of Eq. (I^Tl) in Ref.— involves diagram¬ 
matic perturbation approach and is rather cumbersome. 
However the simplicity of the final result (121]) motivated 
us to re-derive it in a simple and transparent way which 
is presented in the Appendix. 

Let us analyze first Eq. (|2T|) in the inertial interval of 
scales, where according to Fig.jljD, fg = and the terms 
with the viscosities in the denominator may be neglected. 
In this case 


High T > Ta/ 2, two fluid model 



FIG. 4: Color online. Temperature dependence of the ratios: 
Es/Sn - horizontal black line; Rs, Eq. (I28al) . - blue line with 
triangles, Rn, Eq. (|28al) . - green line with diamonds, Rs+n, 
Eq. (I28bl) . - red line with circles. 


Then Eq. (|211) simplihes to the following form, 


ICi{k) ~ 


ap ujg ^ 


ICi{k) 


2/3n7s(fc). 


n -1 


apuJs 

2 pn / k 

a p 


(23) 


(^) 

^ ^max.s ' 


2/3-1 -1 


and it may be analyzed as follows: 

-1 


/Ci(fc) ~ 


1 + 


Pn I’n!' 


ap UJg 


1 -b 


' 1-2 
■ - 


where the viscous cutoff of the superfluid inertial interval 
fcmax.s is given by estimate m- First of all we see that 
the correlation coefficient is governed by the dimension¬ 
less parameter ap/p^ which involves the mutual friction 
coefficient a, as expected. What is less expected, is that 
this parameter, according to Tab.lH depends on the tem¬ 
perature only weakly and is close to unity. Therefore, in 
the inertial interval k <C fcmax.s we have: 


/Ci(fc) ~ 



(24) 


and this expression is very close to unity. In the other 
words, in the inertial interval we expect the full locking 
of the normal and the superfluid velocities for all temper¬ 
atures. This prediction fully agrees with the observations 
in Fig.O 

Consider now case T = 1.07 K, when fcmax.s — 
50 fcniax.n according to the data in Fig.jTlD and esti¬ 
mate (USj). For fcmax.n < k < fcmax.s we have: 

£g{k) ^ £n{k), and > 7 s(fc) > 7 n(^) > . (25) 


Using (dSl), for k ~ fcmax.s we get 


/Ci(A:)~ 


l4/3 


1 -1 


! + ■ 


A /3 


c4/3 

^max, 


« 1, (26) 


We see that the velocities decorrelate in the interval 
^max.n k ^ ^max.s, 3-3 CXpeCted. 

Estimating IC 2 {k) in the regime (l25|) is less simple, be¬ 
cause it requires knowledge of the ratio £g{k)/£n{k) in 
terms of zzn, zz', 7 „(fc) and 7 s(A:). Instead, we can directly 
use Eq. (IA9bl) , which in regime (1251) may be simplified 
(in the (fe, t)-representation) as follows: 

Mn(fe,f) = <-«s(fc,t) . (27) 

p„zz„fc 2 

i.e. Un{k,t) is slaved by Ug{k,t). Equation (07)) immedi¬ 
ately gives /Ci(fc) <C 1, but IC 2 {k) = 1 in full agreement 
with our results in Fig.[3l In particular, this means that 
our simple model of correlations between and Mg, sug¬ 
gested in the Appendix, quantitatively correctly reflects 
the basic physics of this phenomenon. 
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2. Energy dissipation and exchange due to mutual friction 

Strong coupling of the normal and the superfluid ve¬ 
locities suppresses the energy dissipation and the energy 
exchange between the normal and the superfluid com¬ 
ponents caused by the mutual friction (which is propor¬ 
tional to Us — Mn, Eq. (fT^ b Nevertheless, some dissipa¬ 
tion due to the mutual friction is still there. Consider 
the ratio of the total injected energy to the total energy 
dissipated due to the viscosity in the normal and the su¬ 
perfluid components: 


(itn • (itn — Ug)), while the energy dissipation is propor¬ 
tional to the square of this parameter, (|Mn — 


E. Temperature dependence of the effective 
snperflnid viscosity in '^He 

The temperature-dependent effective (Vinen’s) viscos¬ 
ity i''{T) is defined— by the relation between the rate of 
energy-density (per unit mass) flux into turbulent super¬ 
fluid, £, and the vortex-line density, C: 


^s+n — 


Ps^s “f Pn^n 


+ Pnt'n (|w„p) 


(28a) 


a quantity plotted in Fig. S] (red line with circles). Here £„ 
and £s are the inertial range normal and superfluid energy 
fluxes. This ratio exceeds unity by about 10%, meaning 
that ^ 10% of the injected energy is dissipated by the 
mutual friction. As expected, this effect disappears at 
T Ks 1.8 K, when the effective superfluid and normal fluid 
kinematic viscosities are matching (and therefore Us ~ 

lAn). 

The mutual friction has a significantly more important 
influence on the energy exchange between the normal and 
the superfluid components. The energy exchange can 
be quantified by a similar ratio dehned for each fluid 
component. 


Rs = 




Rn = 




(28b) 


shown by a green line with diamonds and a blue line 
with triangles respectively in Fig.lH At the lowest shown 
temperature T = 0.8 K, we have Rn < 0.1 meaning that 
only about 10% of the energy density (per unit mass) 
which is dissipated by the normal fluid component comes 
from the direct energy input. The rest ~ 90% of the 
energy density dissipated by viscosity (at large k) was 
transferred from the superfluid component by the mutual 
friction. This is because for T < 1.8K, we have i^n > v's 
and therefore the normal velocity becomes more damped 
at lower wavenumbers than the superfluid velocity, see 
Fig.Hb- Such an energy transfer by mutual friction from 
the superfluid component to the normal one increases 
Rs (blue line with triangles) above unity. This effect is 
smaller than the one for Rn because at low temperatures 
Ps ^ Pn and the energy per unity volume ps£s + PnCn 
is approximately conserved. As expected, there is no 
energy exchange between the components at T Ri 1.8 K, 
when Vn = t's and Us = Un). At this temperature Rs = 
Rn = Rn+s = 1- Again, as expected for T > 1.8 K, when 
>^n < i^'s (see Tab.|I]) we have Rs > 1, Rn < ^ meaning that 
the energy goes from the less damped normal component 
to the more damped superfluid one. 

To understand why the energy exchange due to the 
mutual friction is larger than the energy dissipation 
by the mutual friction, notice that the energy ex¬ 
change is proportional to the (small) velocity difference, 


£ = z/'(T)(k£)2«z.'(T)(|wsP) . (29) 

According to Eq. (ITSl) . (|wsp) is proportional to the area 
under the plots krt^{k) vs. logfc, shown in Fig. [2] and 
discussed in Sec. Ill Cl These results allow us to determine 
the viscosity v'{T) (analytically and numerically) in the 
entire temperature range from T —>■ 0 up to T —>■ T^. 


1. Low temperature range T < Ta/2 

Consider first the temperature dependence of v'iT) in 
the low temperature range T < T\/2 r; 1.1 K, shown in 
Fig. [5] by the solid blue line. This dependence is found 
in Sec. llVBl in the framework of one-fluid model of grad¬ 
ual eddy-wave crossover Eqs. (ED). As we mentioned, the 
largest (and temperature independent) value of (jwsp) 
(at fixed value of £s) is reached for T < 0.49 K. Ac¬ 
cordingly, the temperature range T < 0.49 K may be 
considered as a zero-temperature limit, at which v'{T) 
reaches its smallest value. The results of the Manchester 
spin-down experiment— are temperature independent as 
well (within the natural scatter of the data). The par¬ 
ticular value of —>• 0) r 0.003k found in these 

experiments is probably accurate up to a numerical fac¬ 
tor (i 3) due to uncertainty in the determination of 
the outer scale of turbulence, taken in Ref— for sim¬ 
plicity as the size of the cube. Our low-temperature, 
one-fluid model (1511) involves one fitting parameter, which 
determines the crossover scale in the blending func¬ 
tion Eq. (|35|) • This parameter affects the resulting value 
of r'(nod(^ 0) chosen such as to meet its ac¬ 

cepted experimental value kp(^ — >■ 0) r 0.003 k. 

At temperatures above 0.5 K, the area under the 
kVil{k) plots in Fig. [2^ become smaller and smaller. This 
is caused by the suppression of the Kelvin wave spectra, 
which is more pronounced at larger temperature, as seen 
in Fig.[D-. The value of (IwsP) decreases with the tem¬ 
perature resulting in a progressive increase in v'{T) as 
shown by the solid blue line in Fig.[5]together with the ex¬ 
perimental (Manchester spin-downi^ and ion-jet—) val¬ 
ues of kp(^)- There is a reasonably good agreement be¬ 
tween the temperature dependence of v'{T) found in the 
framework of the one-fluid model of eddy-wave crossover 
at low temperatures and the experiments. Importantly, 
in the modeling we have used only one phenomenological 
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Low T < Ta/2, one-fluid model High T > Tx/2, two-fluid model 



FIG. 5: Color online. Comparison of the experimental, numerical and analytical results for the temperature dependence of the 
effective kinematic viscosities: Blue triangles - Manchester spin-down experimental^; Green empty circles - Manchester ion-jet 
experiments—; sea-green diamonds with error-bars - Prague counterflow experimental^; cian crosses with error-bars -- Prague 
decay in grid co-flow experiments—; Magenta empty squares - Oregon towed grid experiments2^;Pink right triangles-Oregon 
towered grid experimental. Solid green line - experimental results— for the normal-fluid kinematic viscosity Un = fJ./pn 
(normalized by the normal-fluid density); Dashed green line - He-II kinematic viscosity v = p/p, (normalized by the total 
density) - see also Tab.Hl Thin black dash line - effective viscosity for the random vortex tangle r'(ndi estimated by Eq. (I31II : 
Thick dot-dashed black line - the Vinen-Niemela estimate— of the effective superfluid viscosity, v's, given by Eq. (I32II . Blue 
solid line - v'{T) at T < 1.1 K from numerical solution of Eqs. (ISlbl) for the one-fluid differential model of gradual eddy-wave 
crossover; Red solid line - v' {T) at T > 0.9 K from numerical simulations in Sec. lIII B] of gradually damped two-fluid HVBK 
Eqs. lU in the Sabra shell-model approximation (1331) . 

parameter to fit the zero-T limit of ;Zexp(^)j while the 
temperature dependence of the latter follows from the 
model without any additional fitting. 


2. High temperature range T ^T\l2 

At r > 1 K, Kelvin waves are already fully damped; 
see Fig.[TL. This means that for these temperatures we 
can use the coarse-grained HVBK Eqs. m- Using the 
shell model approximation we find the temperature de¬ 
pendence of r's(T) in the temperature range T > 0.8K 


as shown in Fig. [5] by the solid red line. In the inter¬ 
mediate temperature range 0.9K< T < I.IK this line 
overlaps with the blue solid line, showing the one-fluid 
results. The reason for this overlap is very simple: for 
0.9 K< T < 1.1 K, the Kelvin waves are already damped 
(see Fig.[TL), and the normal-fluid eddies at scales ^ 1 
are still damped. Therefore, in this range both the one- 
fluid model and the coarse-grained model describe the 
physics equally well. Moreover, the effective viscosity v' 
in the one-fluid approximation, suggested in Ref— and 
shown as the black dot-dashed line in Fig. [SI gives the 
same result as the two our approaches. At T < 0.9 K the 
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coarse-grained (blue and black) results deviate below the 
one-fluid prediction which also accounts for the energy 
transfer to the Kelvin waves. This results in a slower 
decrease of v'{T) with temperature, and finally, in the 
zero temperature limit (i.e. below 0.5 K) this predicts 
the plateau which is fully determined by the bottleneck 
energy accumulation at the crossover between the hydro- 
dynamic and the Kelvin wave regimes of the superfluid 
motions. 

For temperatures above 1.1 K, the normal fluid con¬ 
tribution to the two-fluid dynamics becomes important, 
and the two-fluid results in Fig. [5] (red solid line) deviate 
above the value of iz((T) which is determined by the su¬ 
perfluid component alone, cf. Refi^. This is because 
in this temperature interval the normal-fluid viscosity 
t'n(T) (dashed green line) is much higher than its su¬ 
perfluid counterpart Therefore there exists an energy 
flux which is induced by the mutual friction from the less 
damped superfluid component to the normal fluid com¬ 
ponent. This is seen in Fig.lH were i?s > 1 and Rn < 1. 
Thus (|wsP) is suppressed and correspondingly v'{T) de¬ 
viates above i/' (the black dot-dashed line) that does not 
account for the energy exchange. Clearly, at T ~ 1.8K , 
when there is no energy exchange, v' {T) should be equal 
to zz'. Also it is clear that for T > 1.8 K, when the energy 
flows in the opposite direction (from the normal- to the 
super-fluid component), one expects that v'{T) should 
be smaller than v'^. All these expectations are confirmed 
by the results shown in Fig.O 

In Fig.[S] we also show the results for v'{T) of the 
Prague counterflowi^ and co-flow decayi^ experiments 
and those of the Oregon towered grid^^^ experiments. 
These high-T experimental data have a significant scatter 
for reasons discussed in Refs.— “21. 

Taking into account the scatter of the experimental 
data, our computed dependence of v'iT) agrees reason¬ 
ably well with experiments in various flows in the entire 
range of temperatures from T —>■ 0 up to T —>■ T>. 


III. COARSE-GRAINED, TWO-FLUID 
DYNAMICS OF SUPERFLUID TURBULENCE 

This Section concentrates on the high-temperature 
regime, say T > T\/2, when the small-scale motions 
(Kelvin waves) are effectively damped and we can restrict 
ourselves to a coarse-grained description of the super¬ 
fluid dynamics in the continuous-media approximation. 
In most of this temperature range both the normal and 
the superfluid components play important role and the 
two-fluid description is required. 


A. Coarse-grained, two-fluid, gradually-damped 
HVBK equations 

In this subsection we discuss in more details the 
gradually-damped HVBK model, presented by Eqs. (El)- 


1. Simple closure for the large-scale energy dissipation due 
to mutual friction 

Originally in HVBK equations the mutual friction force 
has the form: 

Rns — O: dig X [ctZg X (lAn ^s)] T O: ^ ("^s "^n) • (f^d) 

Here Wg = V x ttg is the superfluid vorticity and the unit 
vector thg = Wg/wg is pointing in the direction of the vor¬ 
ticity. The dimensionless phenomenological parameters 
a and a' describe the dissipative and the reactive mu¬ 
tual friction forces acting on a vortex line as it moves 
with respect to the normal component. 

In our Eqs. ([T]) we used the simplified form (fTB) of the 
mutual friction which accounts for the fact that the vor¬ 
ticity in developed turbulence is usually dominated by 
the smallest eddies in the system, with the Kolmogorov 
viscous scale rj and with the largest characteristic wave- 
vector kr^ ~ 1 / 77 . These eddies have the smallest turnover 
time Trj which is of the order of their decorrelation time. 
On the contrary, the main contribution to the velocity in 
the equation for the dissipation of the fc-eddies with in¬ 
termediate wave-vectors fc, fc A:,j, is dominated by the 
fc'-eddies with fc' ~ fc. Because the turnover time of these 
eddies tj,/ ^ r,,, we can justify the approximation (O 
by averaging the vorticiy in Eqs. (I30|) during the time in¬ 
tervals of interest (r,, <C r ^ Tk'). Thus the vorticity 
Wg may be considered as uncorrelated with the velocities 
Mg, Mn, which are the dynamical variables. More detailed 
analysis22 shows that the approximation (O directly fol¬ 
lows from the Kraichnan’s Direct Interaction Approxima¬ 
tion in the Belinicher-Lvov sweeping-free representation 
for the velocity triple-correlations. 

2. Vinen-Niemela model for the superfluid energy 
dissipation 

The energy dissipation term involving v'^ in Eq. (Hal) 
for the superfluid velocity attempts to take into account 
the existence of quantized vortex lines in an essentially 
classical regime of motion. An early approach to this 
dissipation term— was based on a picture of a random 
vortex tangle moving in a quiescent normal component. 
With the definition (1291) this picture leads to the simple 
equations for the effective viscosity—: 

= (31) 

shown by a black thin dashed line in Fig. [SI One sees 
that this result is much lower than the experimental data. 
Moreover, the picture of a random vortex tangle moving 
in a quiescent normal fluids predicts a higher dissipation 
compared to the realistic situation in which the normal 
and superfluid velocities are almost locked together as 
discussed in Sec. Ill PI The missing physics that needs to 
be considered to resolve this contradiction is that of vor¬ 
tex reconnections. 
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During reconnections, sharp angles necessarily appear 
on the vortex lines, leading to their fast motion. This 
motion is uncorrelated with the motion of other vortex 
lines (except the ones involved in the reconnection) as 
well as with the (relatively slow) motion of the normal- 
fluid component. This leads to large local energy dissi¬ 
pation events due to the mutual friction, which smoothes 
the vortex lines and removes the regions of high curvature 
appearing after the reconnection events. A detailed anal¬ 
ysis of this and related effects led Vinen and Niemela— 
to suggest the effective superfluid (kinematic) viscosity 
which in our notations reads: 

i^i = l3{T)aK, ^(T) = s(|^)'. (32) 

Here, the parameter C 2 relates the vortex line density 
and the mean-square curvatur o^^i^° , and the parameter 
s < 1 accounts for the suppression of effective line den¬ 
sity due to partial polarization of vortex lines and was 
roughly estimated in RefjiS as 0.6. The temperature de¬ 
pendence of r's(T) estimated with the help of Tab. Ill of 
Ref.— is presented in Tab. m The resulting plot of r'((T) 
is shown in Fig.[S] by a thick black dot-dashed line. Be¬ 
sides the clear and definitely relevant physics underlying 
Eq. (I32p . it agrees well with the experiments. That is why 
in our analysis we will include effective damping (13211 in 
our gradually damped HVBK model Eq. (ITal) . 


3. Truncated vs. gradually-damped HVBK models 

A previous model which is referred to as the “trun¬ 
cated HVBK model” of superfluid turbulence was sug¬ 
gested in Ref.—. The idea was to account for the strong 
suppression of Kelvin waves at high temperatures by sim¬ 
ply truncating the HVBK equation for the superfluid at 
a cutoff wavenumber = P/i (at which the normal fluid 
is expected to be well damped by the viscosity), using 
a htting parameter /3 of the order of unity. An obvious 
limitation of this model is the abruptness of the trunca¬ 
tion. An attempt to use this model for the calculation of 
the effective viscosity— shows that the resulting v' [de¬ 
noted in Ref— as z/eff] may vary by a factor of about 
five, when /3 changes by the same factor (see Fig. 4, Top, 
in Ref.—). Moreover, due to the strong temperature de¬ 
pendence of the normal-fluid viscosity z/n, the truncation 
scale depends strongly on temperature, as illustrated in 
our Fig.[T}D. This means that the fitting parameter /3 
should be temperature dependent. If so, the truncated 
HVBK model looses its predictive power. We remark 
however that the experimental values of the effective vis¬ 
cosity presented in Ref.— unaffected by issues with the 
truncated HVBK equations discussed here. 

Unfortunately, this is not the only problem in the anal¬ 
ysis of z/eff in Ref.—. For interpreting the results of the 
numerical simulations of the truncated HVBK model the 
authors of Ref. — use Eq. (|^ for which presumes the 

random vortex tangle. The values of are shown in 


Fig. [5] by a black thin dashed line. As we already noticed 
in Sec. lHI A^ these values are smaller by one order of 
magnitude as compared to the experimental points. The 
physical reason for this discrepancy is very simple. The 
truncated HVBK model ignores the energy dissipation 
in the reconnection events, that, according to our sim¬ 
ulations of gradually-damped HVBK model illustrated 
in Fig.m constitute more than 90% of the total energy 
dissipation in the system. 


B. Two-fluid Sabra shell-model of turbulent '^He 


1. Sabra-shell model equations 


Following— we can present gradually damped HVBK 
Eqs. dH) for isotropic space-homogeneous turbulence as 
the system of two shell model equations for the normal 
and superfluid shell velocities coupled by the friction 
force term. In the dimensionless form it may be written 
as follows: 


d 

dr 

■ d 

-dr 


vnkl, z;" = NL™{z;",}-b—F„-b/"(,33a) 

Pn 




n 


RCk = 




kRck 

LUt 

n 

..S |2 


= NL™{z;^,}-F™ + /^, (33b) 

(33c) 


kRCk ’ 


m = 2 




[Ps-Rs “b Pn-Rn] 




(33d) 


Fm = - v ^), ujI^ = (33e) 


Here NLm{zZm'} is the Sabra nonlinear term, given by 
Eq. (|4|) . The dimensionless shell wave numbers km are 
chosen as a geometric progression km = ^ 02 ™, where m = 
1,2,... M are the shell indices, and the dimensionless 
reference shell wave number is normalized by the inverse 
outer scale of turbulence 1/L (to be specific, in our sim¬ 
ulations we have chosen fco = 1/16). 

The NLm{zzm'} term(j4l) conserves the kinetic energy 
(per unite mass) Ks and ATn (133dl) provided a -I- & -|- c = 0 
(which is our choice: 6 = c = —a/2, a = l). The shell 
energies |z;));®p correspond to the normal- and the super¬ 
fluid energy spectra £n,s{k) as follows: 

|z;));^p = kmSnAkm) , km = . (34) 

Here, the factor km originates from the Jacobian of trans¬ 
formation from dk to dink = (l/k)dk in the integrals for 
the total energy. In particular, the KO-41 spectra with 
^n,s(fc) oc k~^/^ correspond to oc spectra for 

the shell energies. 
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2. Choice of parameters and numerical procedure 


A. Differential approximation for the energy fluxes 
of hydrodynamic and Kelvin wave motions 


A random 5-correlated in time “energy-only” forc- 
ingi^, and /^, was added to the first two shells 
in equations for both the normal- and the super-fluid 
componentsEqs. (1 33 all and (I33bll respectively. Its ampli¬ 
tude was chosen such that the total dimensionless energy 
it's -I- Kn — 1- The standard relation between physical 
velocity fields Un{r), Us{r) and shell velocities 
is a bit involved and will not be displayed here. Note 
only that the dimensionless time r in Eqs. (1331) is nor¬ 
malized by the turnover time of the energy-containing 
eddies. What is more important here is the normaliza¬ 
tion of the viscosities by k and Re^; see Eqs. (1551) . The 
Reynolds number Re^ is a free parameter of the simula¬ 
tions that determines the width of the inertial interval: 
oc Re^'^^. 

The mean effective viscosity in the superfluid subsys¬ 
tem, {v'), is calculated from the mean superfluid energy 
flux, (es), and the mean enstrophy, (wg). According to 
definition (1551) and Eqs. (1551) . we have: 



K 


(gs) 


(35a) 


Here, Wg is given by Eos. (l33eL and the energy fluxes 
through shell m, es.n(fcm), are as follows, 


£s,n(fcm) = fcm Iin[aA53’”(m -I- 1) — cS'g ”(m)](35b) 
Srim) = (35c) 


Eqs. (|551) were solved using the 4th order Runge-Kutta 
method with an exponential time differentiation^!. 

The shell velocities were initiated to have the am- 
plitudes proportional to km and random phases. The 
simulations were carried out for temperatures from T = 
0.8 K to T = 2.15 K using Re^ = 10® and N = 32 shells. 
All other parameters are given in Table ID All observables 
were obtained by averaging over about 500 large eddy 
turnover times. The mean energy fluxes es,n(fcm) were 
calculated by additional averaging over shells 5-^13. 

The results of these simulations are shown in Figs.[T]D, 
[Do, [5113] and E] and were discussed in Sec. Ill A 3l 


IV. LOW-TEMPERATURE, ONE-FLUID 
STATISTICS OF SUPERFLUID TURBULENCE 


In this subsection we discuss the analytic form of the 
energy flux e{k) from small k motions toward the largest 
possible k motions presenting an overview of results for 
£{k) obtained in Refs. 36,61,62,64 required for further 
developments. In Sec. II V AD we begin with the analy¬ 
sis of the expression for e(/c) for the Kelvin wave region, 
M <C I, in terms of its energy spectra Suv/ik) and then, 
in Sec. lIV A~2l we discuss the expression for £{k) for large 
scale hydrodynamic motions in terms of the hydrody¬ 
namic energy spectra. 


1. Small scale motions of Kelvin waves 

It is now recognized that the typical turbulent state 
of a superfluid consists of a complex tangle of quantized 
vortex lines— swept by the velocity field produced by 
the entire tangle according to the Biot-Savart equation—. 
Motions of the superfluid component with characteristic 
scales R i may be considered as motions of individual 
vortex lines, i.e. Kelvin waves. An important step in 
studying Kelvin-wave turbulence was done by Sonin— 
and later by Svistunov— , who found a Hamiltonian 
form of the Biot-Savart Equation for a slightly perturbed 
straight vortex line. The final form of this Hamiltonian, 
found in Ref.— served as a basis for consistent statis¬ 
tical description of Kevin wave turbulence by Lvov and 
Nazarenko— in the framework of standard kinetic equa¬ 
tions for weak wave turbulence— i^. This approach&!!i!! 
led to the spectrum of Kelvin waves E^^ik). Below we 
present a brief overview of these and other pertinent re¬ 
sults. 

a. Zero temperature limit. The total “line-energy 
density” of Kelvin waves (per unit length of the 

vortex line and normalized by the superfluid density) is 
given by the fc-integral of the energy spectrum E.^^(k): 

E’j.w = J E^^{k)dk, E^^(k) = 2uj{k)n{k) .(36a) 

Here n(k) is the wave action, in the classical limit related 
to the occupation numbers Nk as follows: n{k)/h —>■ 
N(k); Lu(k) is the frequency of Kelvin waves. For 
our purposes it is sufficient to use Local-Induction- 
Approximation (LIA)!^ for ui{k)\ 


In Sec. Ill A'Sl we presented an overview of one-fluid de¬ 
scription of superfluid turbulence, based on the differ¬ 
ential approximation for the energy flux in terms of the 
energy spectrum itself. In this section we discuss this 
differential closure procedure in much more details and 
derive second-order ordinary differential equation for the 
superfluid energy spectra £a(T) in the entire range of 
scales, but for T < T^/2. 


uj{k) = k^ , A = In f- ) , (36b) 

47r Vao/ 

A previous model of gradual eddy-wave crossove r^^i®® 
was based on the Kozik-Svistunov (KS) spectrum of 
Kelvin wave turbulence!^ 

KS 

^Kw(*) = '^Ks- -^5 -> KS-spectrum . (37) 
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Here Cj,g ~ 1 is a dimensionless constant and e is the flux 
of (k) in the one-dimensional fc-space. The KS spec¬ 
trum (j37p was obtained in the framework of the kinetic 
equation^ for weakly interacting Kelvin waves under the 
crucial assumptions that the energy transfer in Kelvin 
wave turbulence is a step-by-step cascade, in which only 
Kelvin waves of similar wave numbers effectively interact 
with each other. However in Ref <2^ it was shown that 
this locality assumption is not satisfied. This means that 
KS-spectrum is NOT a solution of the kinetic equations 
and thus physically irrelevant. 

The Kelvin wave turbulence theory was corrected and 
a new local Kelvin wave spectrum was derived by L’vov 
and Nazarenko (LN) in Ref. 

C 

^ fc 5/3 ’ LN-spectrum. (38) 


Here Cln ~ 0.304— and the dimensionless constant 4' is 
given by Eq. (IFbl) . 

This KS vs LN controversy triggered an intensive de¬ 
bate (see e.g. Refs^lvSL)^ which is outside the scope of 
this article. The two predicted exponents, | = 1.4 and 
I ~ 1.67 are very close to each other; indeed vortex- 
filament simulations^^ could not distinguish them (prob¬ 
ably because in this numerical experiment the regime 
of weak turbulence on which the theory is based and 
which requires a small ratio of the amplitude of the 
waves compared to the wavelength, was not the suffi¬ 
ciently realized). Nevertheless, more recent simulations 
by Krstulovio^^, based on the long time integration of 
the Gross-Pitaevskii equations and averaged over an en¬ 
semble of initial conditions (slightly deviating from a 
straight line), support the LN spectrum. The most recent 
vortex-filament simulations by Baggaley and Laurie^! ob¬ 
serve a remarkable agreement with the LN spectrum with 
C™*" ft! 0.308 close to « 0.304 while C"™ « 0.009 
differs from the KS-estimate Cks ~ 1- Based on these re¬ 
sults we will use LN-spectrum ([S]) in further discussions 
of the bottleneck effect. 

b. Differential approximation for the Kelvin-wave 
energy flux. In Ref.—, the LN spectrum of Kelvin 
waves allowed to formulate a differential approxima¬ 
tion for the energy flux, 


£kw if) 


5(ClnAac)3 


dEj^jk) 

dk 


(39a) 


which is an important ingredient of the low-temperature, 
one-fluid differential model. It was constructed by anal¬ 
ogy with Eq. (EH) such as to reproduce the LN spec¬ 
trum (IH together with the thermodynamical equilibrium 
solution £'Kw(fc) =const. The approximation (I39al) plays 
an important role in the discussion of the temperature 
dependence of the effective superfluid viscosity v'{T). 

To analyze the temperature suppression of the Kelvin 
waves energy spectrum, consider the energy balance 
equation 


dcKwik) 


aA 

47r 


— :— nk E/j^\y(/c) 


(39b) 


whose right-hand-side accounts for the dissipation of the 
Kelvin waves in the simplest form suggested by Vinen in 
Ref.—. The approximate solution of Eqs. (IMll . found in 
Ref—, is as follows: 


Ej<iw{k, T) 


Cln a k e. 


1/3 , 


\I/2/3 1 ^ 5/3 

l/£ ^ k ^ k„ 


1 - 




Here cq = eKw(l/^) is the energy influx into system of 
KWs at fc ~ l/£. Notice that all the temperature depen¬ 
dence of Eji^{k,T) is absorbed in fcniax(T) given by: 


^max — ^max (T) 


14^^ 


^/a{T)C^AK 


3/2 


(40b) 


The analytical solution (l40l) is in qualitative agreement 
with the numerical results shown in Fig.[T]j.. 


2. Large-scale hydrodynamie region 

In the hydrodynamic range of scales, the Biot-Savart 
description of superfluid turbulence is too detailed for our 
purposes and we can return to the continuous medium 
approximation, Eqs. o, used in the high temperature 
regime. The difference with Sec. lHIl is that now we 
will first perform the statistical averaging of the velocity 
field, and only then analyze the resulting equations for 
the energy spectra in the one-fluid approximation. This 
different strategy is dictated by a natural requirement 
that Kelvin waves and hydrodynamic eddies have to be 
treated in a similar formal scheme in order to describe 
the intermediate region of scales where one type of mo¬ 
tion continuously turns into the other. As a candidate for 
this scheme we choose a differential closure that allows 
us to express the energy flux as a differential form of the 
energy spectra. For the Kelvin waves this approximation 
was given by Eq. (I39ap and for hydrodynamic eddies it is 
discussed below. 

The simplest approximation for the hydrodynamic en¬ 
ergy flux £HD(fc), based on the Kolmogorov idea of the 
locality of the energy transfer and dimensional reasoning 
goes back to Kovasznay 1947 paper—: 

eHD(fc) ~ C[£^r,ik)f^k^/^ . (41) 

Here C ~ 1 is a dimensionless constant. The basic idea 
of such models is that the nonlinear terms, being of the 
simplest possible form, should preserve the original tur¬ 
bulence scalings and, in particular, predict correctly the 
Kolmogorov cascade. Indeed, in the stationary case and 
in the absence of dissipation the energy flux becomes k- 
independent, £HD{k) => Ehd- Then Eq. (HD) turns into 
the Kolmogorov-Obukhov 5/3-law for SHoik), given by 
Eq. ([H- 

Unfortunately, the simple relation (|4T]) does not de¬ 
scribe the thermodynamic equilibrium (with equiparti- 
tion of energy between states), when the energy flux van¬ 
ishes for ^HD oc This disadvantage is corrected in the 


dk 
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Leith-1967 differential model^, given by Eq. (l7al) . This 
approximation coincides dimensionally with the Kovasz- 
nay model 60, but has a derivative d[Eu-D{k )/that 
guarantees that Enoik) = 0, if Snoik) oc k^. The nu¬ 
merical factor i, suggested in’ , gives the value of the 
Kolmogorov constant Ck 4 i = (24/11)^/^ Ri 1.7 in Eq. dH]) 
that is reasonably close to the experimentally observed 
value. 

A generic hydrodynamic spectrum with a constant en¬ 
ergy flux was found in’ as a solution to the equation 
Enoik) = Ehd = const: 


2/3 

Suoik) = Ck41-^ 


1-f 


f—1 

keq' 


1 2/3 


(42) 


in which kgq is as yet a free parameter describing the 
crossover between the low-fc K041 spectrum (jb)) and the 
thermalized part of the spectrum, Snuik) oc with 
equipartition of energy at large k. 

Notice, that Eq. does not account for the energy 
dissipation due to the mutual friction and the viscosity. 
We will do this later, introducing dissipation terms in the 
energy balance equation and numerically solving them. 


B. One-fluid differential model of gradual 
eddy-wave crossover 

1. “Line-” and “volume-” energy densities, spectra and 
fluxes 

Our goal here is to formulate a model which will al¬ 
low to describe in a unified form the hydrodynamic en¬ 
ergy flux e:HD(fc) at small k and the corresponding ob¬ 
jects for the Kelvin waves. However we cannot do it 
straightforwardly using the equations for Cnwik) and the 
spectrum Ey^^ik). The reason is simple: the hydrody¬ 
namic and Kevin wave objects have different physical 
meaning and different dimensions. Indeed, the hydrody¬ 
namic motions fill the three-dimensional space (volume), 
their energy density £1 hd per unit mass has a dimension 
[^HD]=cm^/sec^. Accordingly, the dimensions of energy 
spectrum £-ay,{k) and the energy flux En-D(k) are as fol¬ 
lows: 

[fHD(fc)] = cm^/sec^ , [£:HD(fe)] = cm^/sec^ . (43a) 

On the other hand, Kelvin waves propagate along one¬ 
dimensional lines - vortex filaments. Therefore the en¬ 
ergy density of Kelvin waves on individual vortex fila¬ 
ment Ekw is normalized by unit vortex length and (for 
the sake of convenience) by superfluid density. Therefore 
its dimension is [EKw]=cm^/sec^. Then the dimensions 
of the corresponding energy spectrum and the energy flux 
are: 

[Ey^{k)] = cm®/sec^ , [eKw(fc)] = cm^/sec^ . (43b) 

Different normalization of the same objects dictates 
the relation between them in a statistically homogeneous 


and isotropic vortex tangle with the line density L = l/£^ 


4w(fc) 


-^Kw(fc) 

£2 


^KW 



(44) 


2. Energy balance equation 

Consider a general form ([8]) of the continuity equa¬ 
tion for the energy density £s(k, t) of the isotropic space- 
homogeneous turbulence of the superfluid component 
which accounts for the energy dissipation with the help 
of the Vinen-Niemella viscosity. The remaining phys¬ 
ical problem here is how to describe the energy density 
£s{k, t), the energy flux over scales Es{k, t), and the damp¬ 
ing term 'D{k,t) in the entire range of wave vectors fc, 
including the intervortex scales k ^ IJL A step toward 
this direction was suggested in Ref. [ll| in the form of 
the “Eddy-wave model” in which the superfluid motions 
with scales R ^ £ are considered as a superposition of 
two coexisting and interacting types of motion: random 
eddies and Kelvin waves. In some sense the problem 
here is similar to the description of the mechanics of the 
matter at intermediate range of scales, where it behaves 
like particles and waves simultaneously. In quantum me¬ 
chanics it was suggested to formulate explicitly the basic 
equation of motion (the Schrodinger equation) and to 
compare its prediction with observations. We are not 
so ambitious, our goal is to discuss below a set of un¬ 
controlled approximations, based currently only on our 
physical intuition, which leads to an explicit set of model 
equations having a predictive power. As we will see be¬ 
low, our model predicts that in the range R ^ i, due 
to the bottleneck energy accumulation, the energy dis¬ 
tribution between scales is close to the energy equipar¬ 
tition, like in the thermodynamic equilibrium. It is well 
known that in thermodynamic equilibrium the statistics 
is universal and independent of the details of interac¬ 
tion. Therefore we hope that many details of the vortex 
dynamics (including the vortex reconnections), that are 
ignored in our model, do not affect the model results: we 
believe that these results are closer to the reality than 
the model itself. 


3. Gradual model for the energy spectra of superfluid 
component 

The basic physical idea is to approximate the total tur¬ 
bulent superfluid energy density £s{k,t) as a sum of the 
hydrodynamic energy spectrum £^y,{k,t) and the energy 
spectrum of the Kelvin waves £Kw{k,t), with the energy 
distribution between the components depending only on 
the dimensionless blending function g{k£) of the dimen¬ 
sionless wave-number ki\ 

£s{k,t) = £^y,{k,t) -\- Eyvfikfl), 

^nnikfl) = g(ki)£s{k,t), (45) 

£Kw{k,t) = [1 - g{k£)]£s{k,t) . 
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In order to find a qualitative form of the blending func¬ 
tion gikt) we follow Refii. Consider a system of locally 
near-parallel vortex lines (in the vicinity of some point 
ro), separated by the mean distance Denote the indi¬ 
vidual vortex lines by an index j. Notice that in principle 
the same vortex line can go far away and come close to vq 
several times. To avoid this problem one should assign 
the same vortex line a different index j if it leaves (or en¬ 
ters) the ball of radius i^/A centered at Tq. Each vortex 
line produces a superfiuid velocity field Mj (r), which can 
be found by the Biot-Savart Law. 

The total superfluid kinetic energy density (per unit 
mass) £s = | Si i (“i 'be divided into two 
parts, Es = Sis + S 2 s, where 


however the fluxes and e-nw are not equal to the fluxes 
£r|ij 3 , Eq. (ITal) . and etkw in isolated hydrodynamic and 
Kelvin wave systems. Equation (l39aD for eKw(A) in the 
volume normalization takes the form (I7b|) . 

The fluxes and ?kw contain additional cross¬ 
contributions (k) and that originate from the 

interaction of two types of motion, hydrodynamic and 
Kelvin waves: 

KW ^ HD 

^HD ^HD ^HD '! ^KW ^KW ^KW * (49lDj 

We modeled the cross-terms in the linear approximation 
with respect to the energies (i.e. the Hydrodynamic en¬ 
ergy affecting the Kelvin waves flux and vice versa): 


<«) 

3 

i<j 

The same subdivision can also be made for the en¬ 
ergy spectrum in the (one-dimensional) fc-space, Ss{k) = 
Sis{k) -\-S 2 sik), with two terms, that may be found via k- 
Eourier components of the superfiuid velocity fields t)|(fc) 
similar to Eq. (H51) . Now the idea is as follows: the energy 
£is{k) is defined by the form of the individual vortex lines 
that is determined by the Kelvin waves, while the energy 
S 2 s{k) depends on correlations in the form of different 
vortices, that produce collective, hydrodynamic type of 
motions. Therefore £is{k) may be associated with the 
Kelvin waves energy, fis(fc) £nyj{k), while £ 2 &{k) has 
to be associated with the superfiuid hydrodynamic en¬ 
ergy, £ 2 s{k) => S^^ik). This allows us to conclude that 


^HD (k) = d [£^^{k)/kl^/dk'^ , (49c) 

4w(fc) = -t^kwI^kw} d[si^{k)/e]/de , 

where ken is the wave number, at which g{ken(-) = \- 
The differential form of these contributions follows from a 
physical hypothesis that these terms should disappear (or 
become negligibly small) when the influencing subsystem 
is in thermodynamical equilibrium, i.e. when oc 
and ^Kw oc kP =const. Eunctionals of the correspond¬ 
ing energies, may be modeled by dimensional 

reasoning, in the same way as Eqs. (I7al) and (I39a|) were 
formulated for the fluxes. The resulting equations for 
JnD and J^kw may be written in the form: 

J-hd{£L} = (49d) 

-^Kwl^Kw} = C'Kw(fc^) • 

Here 


g{kl)=[l+£is{k)/£ 2 s{k)\ \ (47) 

The rest are technical details presented in Ref. [a where 
it was concluded that in practical calculations it is rea¬ 
sonable to use an analytical form g{k£) of the blending 
function 


g{ke) 

goiki) 


go[ci ln(A -|- 7.5) ki] , 



(M)2 exp(A:£) 
47r(l -I- k£) 


(48) 


where ci ~ 0.32 is the fitting parameter, chosen inii. 


4- Gradual model for the superfluids energy flux £s{k,t) 

Modeling the total superfiuid energy flux over scales, 
£s{k, t), is less straightforward than the model (1451) for the 
energy itself, Ss{k,t). We first assume that £s{k,t) may 
be presented as the sum of the fluxes over hydrodynamic 
and Kelvin wave components, 

Es = e^D + ?KW , (49a) 


Chd = -1/8 , C'kw(H) = -5(H)®/7A® . (49e) 


as explained in Ref.—. The resulting model for the total 
energy flux e(fc) follows from Eqs. (j49|): 


£s{k) = -{g\/fc“5(fc£)fs(fc) 

, 3 4/2 (kefkUi-g{ki)?Ssik)^ 


-{ 

dk\ 


Ss{k) 


(C',hA/c)3 
g{k£) , l-g{k£) 


L A:2 


i-2 

'/en 


} (50) 


Only with the choice (I49el) the resulting Eq. (ISOl) for £s{k) 
vanishes in thermodynamic equilibrium (with £s{k) oc 
k'^ in the hydrodynamic regime, k < ken and with 
Ss{k) =const. in the Kelvin waves regime, k > ken,) 
as required. 


5. Dimensionless form of the gradual one-fluid model 

The resulting Eqs. ([5]), (H51) . (H51) and (15(11) represent our 
eddy-wave model of superfiuid turbulence in the one-fluid 
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approximation, which neglects motions of the normal- 
fluid component, assuming = 0. Now we introduce 
dimensionless variables: 

X = M, e{x) = —— , e(x) = —^ (51a) 

and express Eq. (1501) in a dimensionless form which is 
convenient for numerical analysis: 


d rl 


dx 


;\/ 9 {x)e[x) + 


3 [^'a:ena;^(l - g[x)e{x))\ ‘ 




dx 


{e(a 


9{x) , l-9{x) 


]} 


a^ojT-gix) + x'^P{T)'^e{x), 

Xen = keni — 6.64/ ln(A -|- 7.5) . 


(51b) 


(51c) 


In the dimensionless form there is a constraint on the 
energy flux: 

noc 

2 / x^g{x)e{x)dx = 1 , (51d) 

Jo 

following from the assumption that the vorticity is dom¬ 
inated by the scales of order Xjl, given by Eq. (31b) in 
Ref.—. The function defined by Eq. (|5bl) in dimension¬ 
less variables should be found self-consistently by enforc¬ 
ing the following condition, 

4' = —/ [I - g{x)]e{x)dx , g{xi/ 2 ) = - ■ (51e) 

•lxx/2 ^ 


6. Numerical procedure 


We solved the integro-differential Eq. (I51b|) numeri¬ 
cally starting from the large x region. 

To formulate two boundary conditions at large x, we 
use an analytical form of the Kelvin waves spectrum. In 
dimensionless form Eq. (I40al) reads: 


Ckw)^ t Xniax) 


^max 


C 

^2/3 j.5/3 


14^^ 

(aC',,A2)3/4 ■ 



(52a) 

(52b) 


Now we can take as the boundary conditions the values 
of eKw{x) at two points, e(xmax - xi) and e(a;niax - X 2 ) 
with some appropriate values of xi and X 2 (say 5 and 
10) for very small a. The results of these simulations are 
shown in Eigs.H^, [2^, and [5] and were discussed in Sec. HI 


V. SUMMARY AND DISCUSSION 

In this paper, we have generalized the zero tempera¬ 
ture theoryiii2& of the energy and the vorticity spectra 


in superfluid turbulence to non-zero temperatures up to 
T —>■ Ta, accounting for the effect of the mutual friction 
and motion of the normal fluid component. In particular 
we describe the influence of the temperature on the bot¬ 
tleneck energy accumulation near the inter-vortex scales. 

• The gradually damped HVBK Eqs. ([T]) include the 
Vinen-Niemela superfluid viscosity (I32|) with a fitting pa¬ 
rameter s ~ 0.6 which was chosen in their paper—. Be¬ 
sides this, our Sabra-model Eqs. (1^ which is based on 
Eqs. © and is used in the T > Tc/2 range, has no addi¬ 
tional fitting parameters. 

• The differential one-fluid model of superfiuid tur¬ 
bulence Eqs. (I51bl) . used in the T < Tc/2 range, has 
only one fitting parameter ci ~ 0.32, entering into the 
blending function (l4^ and chosen in Ref.—. Besides this 
the model has no additional fitting parameters. Thus in 
the entire approach we used only two fitting parameters 
which were chosen in previous papers. 

• We have shown that for T < 0.5 K Kelvin waves are 

excited in the range of scales from ~ 1 up to some tem¬ 
perature dependent cutoff fcmax (140bl) : see Eig.[TJa,. For 
k£ > 20 Kelvin waves have the LN-energy spectrum ([5]) 
with a constant energy flux, oc while in 

the crossover region (about one decade around ~ 1) 
there exists a flux-less spectrum £^^(k) Riconst corre¬ 
sponding to the thermodynamic equilibrium with the en¬ 
ergy equipartition between Kelvin waves with different k. 
In this temperature range the effective superfiuid viscos¬ 
ity may be considered as temperature independent and 
equal to its zero-temperature limit, v'{T) Ri !z'(0). Also, 
a minor amount of the normal-fluid component may be 
completely ignored. 

• When T exceeds ~ 0.5 K, the constant energy flux 
range of Kelvin waves disappears and the flux-less range 
(with £’kw(^) ~const) begins to shrink; see Fig.[T^. This 
leads to the temperature suppression of the bottleneck 
energy accumulation. As a result, the superfiuid square 
vorticity, (Iwgp), decreases (see Fig. [2^), leading to in¬ 
crease in v'{T)^ in accord with the experimental obser¬ 
vations; see Fig.jSJ Up to T ~ 1.1 K, a small amount 
(below 1%) of the normal-fluid component may be con¬ 
sidered as being at rest, at least in the region M ^ 1, 
which determines the leading contribution to (jo’sP). As 
a result, both our models (the one-fluid gradual model of 
the bottleneck crossover, which account for the presently 
negligible energy flux by Kelvin waves, and the gradu¬ 
ally damped HVBK model, that ignores Kelvin waves, 
but accounts for the presently negligible normal-fluid 
motions) are valid for 0.5K ;< T < IK. The temper¬ 
ature dependence of y'{T) predicted by these models 
(red and blue solid lines in Fig. [5]) practically coincide 
for 0.5K<T< 1.1 K. 

• In the high temperature regime, T > T\/2 ~ I.IK, 
the normal fluid component begins to play some role in 
the temperature dependence of y'{T). In spite of the 
almost full super- and normal-fluid velocity locking (see 
Fig. [3]) there is a significant energy exchange between the 
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components (see Fig. 11, caused by the mutual friction 
and the velocity decorrelation near k£ ^ 1. This physical 
effect is important by itself, although it leads only to a 
small, (but visible) deviations of the resulting tempera¬ 
ture dependence of i^'iT) (solid red line in Fig.H from 
the Vinen-Niemella model^ of i^'(r) = Eq. (|^ . 

shown by black dot-dashed line. 

• Since there is no detailed information on the su¬ 
perfluid and normal fluid energy spectra, especially for 
large k and low temperatures (see e.g. review^), we 
compare our results with the experiments for the tem¬ 
perature dependence of the effective kinematic viscosity 
v' {T). The latter was measured in the temperature range 
from 0.08 K to 2.15 K by the Manchester spin-dowi*i^ and 
ion-jet^S experiments, as well as the Oregon towed-grid^S 
and the Prague counter-flowi^ experiments, all shown in 
Fig. [51 Our computed temperature dependence of the 
effective viscosity v'{T) agrees qualitatively with the ex¬ 
perimental data in the entire temperature range: from 
T —>• 0 up to r T\. We consider this agreement as 
a strong evidence that our low-temperature, one-fluid 
differential model and high temperature coarse-grained 
gradually damped two-fluid HVBK model capture the 
relevant basic physics of the turbulent behavior of ^He. 

•The models considered in this paper are intended 
for systems whose anisotropy effects are not substantial. 
Strong external rotation may change the behaviour by 
enforcing a strong polarisation of the superfluid vortex 
bundles leading to the suppression of reconnections. In 
turn, suppressed reconnections result in an enhanced bot¬ 
tleneck accumulation of the turbulent spectrum near the 
crossover scale £ and, as a result, in a decrease of the ef¬ 
fective viscosity v'. We leave the study of such an effect 
of strong polarisation on the bottleneck phenomenon to 
future. 

• In this paper, we have ignored the effect of mu¬ 
tual friction on the small k region of scales in the low 
temperature regime, when the normal component is rare 
and motionless due to a very large kinematic viscosity 
(i^n > 38 K for r ^ IK). This is justified when the range 
of scales greater that £ is not very wide, as in all existing 
^He experiments. Theoretically, the mutual friction effect 
grows as k is decreased and, if the \ow-k range is wide, 
the spectrum would inevitably reach a friction-dominated 
scaling regime with a power-law exponent equal —3, see 
Ref.—. Such a regime, which is even more natural in ^He 
turbulence, may lead to a vortex tangle decay law with 
L ^ 1/t. We leave the study of this dissipative regime 
for the future. 
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Appendix A: Simple model of cross-velocity 
correlations in superfluids 

Our goal here is to suggest a relatively simple and phys¬ 
ically transparent model of the cross-correlation function 
of the normal and superfluid velocities, that leads to 
Eq. ([211) in the simplest case of homogeneous, isotropic 
turbulence of incompressible turbulent motions of ^He. 

To start, we recall some definitions and relationships, 
required for our derivation, which are well-known in sta¬ 
tistical physics. The first one is Fourier transform in the 
following normalization: 

/ dk 

Un,sir,t)exp{-ik ■ r)] . (AI) 

Next we define simultaneous correlations and cross¬ 
correlations in the fc-representation, [proportional to 
S{k~ k') in homogeneous case]: 

(u„(fc,t)-<(fc',t)) = (27r)3Gn„(fe)J(fe-fe'),(A2a) 
(us(fc,t)-<(fc',£)) = (27r)3G,,(fc)(5(fc-fc'), (A2b) 
(u„(fc,f)-<(fc',t)) = (27r)3Gns(fe)(5(fc-fc') ■(A2 c) 

It is known that their fc-integration produces one-point 
correlations: 


/ dk 

j^Gnnik) = {\un{r,t)\^) , (A3a) 

/ dk 

= (|Ms(r,t)p), (A3b) 

/ dk 

^ 2 ;^Gns(fc) = {Un{r,t) ■ Us{r,t)) . (A3c) 

In isotropic case, each of three correlations G...(fc) is in¬ 
dependent of the direction of k: G,,,{k) = G,,,{k) and 
J .. .dk = 4 tt J .. .k"^ dk. Together with Eqs. (1181) and 
(IA3I) this gives: 


£n{k) = ^Gn„(fc), fs(fc) = I^Gss(fc), 

£^,{k) = ^Gns(.k) . (A4) 

To begin with the derivation of Eq. (I^Tl) . we simplify 
Eqs. dH) for the superfluid and the normal velocities, 
Vs(k,t) and Un(fe,£), by modeling the nonlinear terms 
in the spirit of the Langevin approach, i.e. replacing 
them by a sum of respective damping terms 7 sUs(/i:,t) 
or 7nMn(fc,t) and random, delta-correlated in time force 
terms fs{k,t) or fn{k,t) with Gaussian statistics and 
zero cross-correlations: 

{Uk,t)-f:{k',t')) = {2nrSik-k')S{t-t')fi{k), 
{fn{k,t)-f:{k',t')) = {2^f5{k-k')5{t-t')fl^{k), 
(/s(fc,t)./:(fc',0) = 0 . (A5) 
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In the fc-representation, the resulting equations read: 


+ TsVs{k,t) = auJsVn{k,t) + fs{k,t), (A6a) 

dVn{k,t) YnVnik,t) = —WsVs(fc, t) +/n(A;, t), (A6b) 
at Pn 

hn = 7n +, Fs = 7s + aws + ■ (A6c) 

Pn 


Multiplying Eqs. (lABal) and (IA7bl) by Vg, and Vn, respec¬ 
tively and averaging, one gets equations for velocity cor¬ 
relations Gnn, Gss and cross-correlation Gsn, defined by 
Eqs. (IA3I) : 


-_a 

.m 

rd 

.m 


9 

1^ 


-2r, 

2r„ 

Fs -|- Fn 


Gs: 

G„. 

Gs, 


2aujgGgn + 2Re[$ss] , 

(A7a) 

2—^WsGsn + 2Re[$nn] , 

(A7b) 

Pn 

- \aps^ ^ 1 

Ws - Gnn + CtsGss + 

(A7c) 

1 Pn J 

Re[$sn + $ns] ■ 


These equations involve yet unknown cross-correlations 
of the velocities and the forces, defined similarly to 
Eqs. (IA2I) : 

(/„(fc,t).<(fc',t)) = (27r)3$„„(fc)J(fc-fc'), (A8a) 
{fs{k,t) ■ v*{k',t)) = (27r)^d>ss(fc) - fc'), (A8b) 

(/„(fe,t)-<(fc',t)) = (27r)3$,,(fc)J(fc-fc'), (A8 c) 

(/s(fe,t) •<(&', t)) = (27r)3$,„(fc)J(fc-fc') . (A8d) 

To find these correlations, we rewrite Eqs. (IA6I) in Fourier 
w-representation: 

[iw-I- Fs]us(fc, w) = asCdsVn{k,uj) + fg{k,uj), (A9a) 

[itu + rn]vn{k,uj) = — WsS's(A;,a;)-b/n(fe,u;) , (A9b) 


were u,,. and f . denote Fourier transforms of the cor¬ 
responding functions. The solution of linear Eqs. (IA9I1 
reads: 


Vg = -[{iuj + rn)fs + C(gUJgfn]/A, (AlOa) 

Vn = -[(iw-b Fs)^-b OsWs/sj/A , (AlOb) 


A = 


2 , ^ Ps -2 -T' -T' 

OJ H- UJ — 1 si n 


iw(Fs -b Fn) , (AlOc) 


where for brevity we skipped the arguments {k,uj) in all 
functions. Multiplying the two Eqs. (lAlOl) by /„ and 
fg, respectively and averaging, we get four equations 
for (cross)-correlations $ (fc, w) in ^-representations via 
two correlations fgg{k) and /nn(fc) of random forces, 
Eqs. (IA5I) . By integration of the results over w one may, 
in principle, get the simultaneous cross-correlation func¬ 
tions $nn(fc), ‘f’ss(fe), ‘hsn(fc) and $ns(fc), cxprcsscd via 
fgg{k) and /^„(fc). However, it suffices to realize that 
^sn{k) = ^ns{k) = 0. This instantly simplifies the sta¬ 
tionary (d/dt = 0) Eq. (IA7cl) to 


Ts + Tn 

Gsn — 

^nn “T UsLxgs 



Pn J 


(All) 


and we get a relationship between Gnn(fc), Ggg{k) and 
Gsn(fe), that is equivalent [with account of Eqs. (EH)] to 
Eq.®. 
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